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Abstract 

The production of energetic neutrons in ep collisions has been studied with the 
ZEUS detector at HERA. The neutron energy and distributions were mea- 
sured with a forward neutron calorimeter and tracker in a 40pb~^ sample of 
inclusive deep inelastic scattering (DIS) data and a 6 pb"^ sample of photopro- 
duction data. The neutron yield in photoproduction is suppressed relative to 
DIS for the lower neutron energies and the neutrons have a steeper distribu- 
tion, consistent with the expectation from absorption models. The distributions 
are compared to HERA measurements of leading protons. The neutron energy 
and transverse-momentum distributions in DIS are compared to Monte Carlo 
simulations and to the predictions of particle exchange models. Models of pion 
exchange incorporating absorption and additional secondary meson exchanges 
give a good description of the data. 
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1 Introduction 



In ep scattering at HERA, a significant fraction of events contains a low-transverse- 
momentum baryon carrying a large fraction of the incoming proton energy [1-5]. Al- 
though the production mechanism of these leading baryons is not completely understood, 
exchange models (Fig. [1]) give a reasonable description of the data. In this picture, the 
incoming proton emits a virtual particle which scatters on the photon emitted from the 
beam electron. The outgoing baryon, of energy Eb, carries a fraction xl = Es/Ep of 
the beam energy, while the exchanged particle participates in the process with energy 

(1 - XL)Ep. 

In particular, one-pion exchange is a significant contributor to leading neutron production 
for large xl [1,2]. For such a process the cross section for the semi-inclusive reaction 
7*p —>■ Xn factorizes into two terms (Regge factorization [6]): 

-^^^^ = fn/p{xL, t)a^*^{{l - xl)W , Q ), 

where is the virtuality of the exchanged photon, W is the center-of-mass energy of the 
virtual photon-proton system and t is the square of the four-momentum of the exchanged 
pion. In terms of the measured quantities Xl and transverse momentum pt, the pion 
virtuality is: 



t ~ - 



pI {l-XL){ml-mlxL] 



Xl Xl 

The flux of virtual pions emitted by the proton is represented by /^r/p and ct^.tt is the 
cross section of the virtual-photon and virtual-pion interaction at center-of-mass energy 
^/l — Xl W. If the 7*7r cross section is independent of t, the pr distribution of produced 
neutrons is completely determined by the flux factor. 

Many parameterizations of the pion flux have been suggested in the literature [6-11]. 
They have the general form: 

L/p{xL,t) OC ~\ (l-XLT'-'^F^XL^t). 

The power a{t) and the form factor F{xL,t) are model dependent with parameters that 
can be extracted from hadron-hadron scattering data. 

Comparisons between cross sections for the production of particles in the fragmentation 
region of a target nucleon provide tests of the concepts of vertex factorization and limiting 
fragmentation [12]. The hypothesis of limiting fragmentation states that, in the high- 
energy limit, the production of particles in the proton target-fragmentation region is 
independent of the nature of the incident projectile. For leading neutron production in 
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ep scattering, where the projectile is the exchanged virtual photon, this implies that the 
dependence of the cross section on the lepton variables {W, Q"^) should be independent 
of the baryon variables {xL,t). For such vertex factorization, the cross section can be 
written as 



where g and G are arbitrary functions. The Regge factorization introduced earlier violates 
this vertex factorization because of the dependence of a^*-,^ on xl in addition to 
and Q^. Factorization tests involve comparing semi-inclusive rates, normalized to their 
respective total cross sections, to study whether particle production from a given target 
is independent of the lepton variables. 

In exchange models, neutron absorption can occur through rescattering [13-16]. In a 
geometrical picture [14] , if the size of the n-vr system is small compared to the size of the 
photon, the neutron can also scatter on the photon. The neutron may migrate to lower 
Xl and higher such that it is outside of the detector acceptance. The rescattering 
can also transform the neutron into a charged baryon which may also escape detection. 
Since the size of the virtual photon is inversely related to Q^, more neutron rescattering 
would be expected for photoproduction (Q^ ~ 0) than for deep inelastic scattering (DIS, 
^ 1 GeV^). A previous study [2] showed a mild violation of vertex factorization with 
the expected increase of rate when going from photoproduction to DIS. Similar effects 
have also been seen for leading protons [3]. The size of the n-vr system is inversely 
proportional to the neutron pt, so rescattering removes neutrons with large pt- Thus 
rescattering results in a depletion of high p^ neutrons in photoproduction relative to DIS: 
a violation of vertex factorization. Pion-exchange models [6-11] incorporate a variation 
of the mean size of the n-vr system as a function of x^. This results in an dependence 
of the absorption, again a violation of vertex factorization. 

Absorption is a key ingredient in calculations of gap-survival probability in pp inter- 
actions at the LHC, critical in interpreting hard diffractive processes, including central 
exclusive Higgs production. The most recent absorption model calculations [15,16], based 
on multi-Pomeron exchanges, gave a good description of previous leading-neutron results 
on absorption [2]. 

This paper presents measurements of the xl and p"^ distributions of leading neutrons 
coming from samples of DIS and photoproduction (7p) processes, with more than seven 
times higher statistics and smaller systematic uncertainties than the previous ZEUS pub- 
lication [2]. The Xl and distributions in DIS and photoproduction are compared as 
a test of vertex factorization. The neutron measurements are compared to similar mea- 
surements of leading protons at HERA. The data are also compared to the predictions of 
several Monte Carlo (MC) models. The neutron distributions in DIS are compared to 
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several pion-exchange models with various choices of their parameters. Finally, the xl and 
Pt distributions in photoproduction and DIS are compared to models incorporating pion 
exchange and rescattering, and a model that also includes secondary meson exchanges. 

2 Detectors 

A detailed description of the ZEUS detector can be found elsewhere [17]. A brief outline 
of the components that are most relevant for this analysis is given below. 

The high-resolution uranium-scintillator calorimeter (CAL) [18] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = O.IS/a/E for electrons 
and a{E)/E = 0.35 /\/E for hadrons {E in GeV). The EMC sections were used to detect 
scattered positrons in DIS events and the RCAL was used to trigger on the dissociated 
photon in photoproduction events. 

Bremsstrahlung, ep — ^ ejp, and the photoproduction of hadrons, ep eX, are tagged 
using the luminosity detectors [19]. The bremsstrahlung photons are measured with a 
lead-scintillator calorimeter located at Z = —107 nfl from the interaction point in the 
positron-beam direction. The positron tagger was a similar calorimeter at Z = —35 m 
from the interaction point with an energy resolution of a(E)/E = 0.19/V^ {E in GeV). 
It was used to measure positrons scattered at very small angles in an energy range of 
5-20 GeV. 

The forward neutron calorimeter (FNC) [20] was installed in the HERA tunnel at ^ = 
degrees and at Z = 106 m from the interaction point in the proton-beam direction, as 
depicted in Fig. [21 It was used for the 1995-2000 data taking. The FNC was a lead- 
scintillator calorimeter with an energy resolution for hadrons measured in a test beam 
to be a{E)/E = 0.70/y/E (E in GeV). The calorimeter was segmented vertically into 
14 towers as shown in Fig. [31 Three planes of veto counters were located in front of the 
FNC to reject events in which a particle showered in dead material along the beamline 
upstream of the FNC. 

In 1998 a forward neutron tracker (FNT) was installed in the FNC at a depth of one 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
center of HERA. The coordinate origin is at the nominal interaction point. 
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interaction length. It was a scintillator hodoscope designed to measure the position of 
neutron showers. Each scintillator finger was 16.8 cm long, 1.2 cm wide and 0.5 cm 
deep; 17 were used for X position reconstruction and 15 for Y. Figure [3] shows the 
position of the FNT hodoscope in the FNC relative to the incoming neutron beam. The 
irregular outlined area indicates the geometric acceptance defined by magnet apertures. 
This limited detection to neutrons with production angles less than 0.75 mrad, allowing 
transverse momenta in the range pt < En ^max = 0.69 xl GeV. The resulting kinematic 
regions in and t are shown in Fig. HI 

Scans by a ^'^Co radioactive source and data from frequent proton beam-gas runs were 
used to calibrate and monitor both detectors. The relative calibration between FNC 
towers was adjusted using position information from the FNT. The energy scale of the 
FNC was determined with a systematic uncertainty of ±2% from fits to the endpoint of 
the neutron energy spectrum near 920 GeV. The minimum-ionizing-particle (mip) scale in 
the veto counters was determined by selecting electromagnetic showers in the FNC, a large 
fraction of which were converted photons which deposited a 2-mip signal in the counters. 
The position resolution of neutron showers in the FNT of ±0.23 cm was measured by 
placing an adjustable collimator in front of the outermost veto counter of the FNC during 
special test and calibration runs. 

3 Data selection and analysis 

The data for this analysis were collected in 2000 when HERA collided 27.5 GeV positrons 
with 920 GeV protons, giving an ep center-of-mass energy ^/s = 318 GeV. Separate trig- 
gers were used to collect DIS and photoproduction events with leading neutrons. 

3.1 Data selection 

The DIS events were collected using a trigger that required the detection of the scattered 
positron in the CAL. In the offline analysis, the scattered positron was required to have 
energy E'^ > 10 GeV and to be at least 3 cm from the inner edge of the beam-pipe hole 
in the RCAL. The quantity E — Pz = '^iEi{l — cos^j), with the sum running over all 
calorimeter cells, was required to be in the range 35 < — < 65 GeV; the lower cut 
reduced photoproduction background with a misidentified positron. These cuts resulted 
in a clean sample of DIS events in the kinematic range Q2 > 2 GeV^ with a mean photon 
virtuality of (Q^) — 13 GeV^. For further studies, the variable was reconstructed using 
the double angle (DA) method [21]. This method requires a certain amount of hadronic 
activity in the CAL in order to measure the angle of the hadronic system. To ensure 
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this, an additional requirement yjB > 0.02 was imposed for those measurements requiring 
reconstruction. Here ?/jb is the inelasticity y ^ [W"^ + Q'^)/s reconstructed using 
the Jacquet-Blondel (JB) method [22]. The integrated luminosity of the DIS sample was 
approximately 40pb~^. 

The photoproduction events were collected during the last part of the 2000 running period 
using a trigger that required at least 5 GeV in the positron tagger in coincidence with at 
least 464 MeV in the RCAL EMC [23]. The acceptance of the positron tagger limited 
the photon virtuality to < 0.02 GeV^, with a mean of approximately (Q^) ~ 
4 X 10"^ GeV^ Offline, the total energy per event deposited in the photon tagger was 
required to be less than 1 GeV in order to reject overlapping bremsstrahlung events. The 
integrated luminosity of the photoproduction sample was approximately 6pb~^. 

The DIS and photoproduction triggers required at least 180 GeV of energy to be deposited 
in the FNC. Good FNC neutron candidates were required offline to satisfy the following 
conditions: 

• the FNC tower with maximum energy was one of the four towers covered by the FNT 
as depicted in Fig. [31 to reject protons with xl < 1 which were deflected into the top 
towers of the FNC by the vertical bending magnets shown in Fig. [2l 

• the veto counter had a signal of less than one mip, to reject showers which started in 
dead material upstream of the FNC; to minimize effects of backsplash from hadronic 
showers, only the veto counter farthest from the FNC was used; 

• no signal in the veto counter consistent with a shower from a previous bunch crossing, 
to reject pile-up energy deposits; 

• the timing information from the FNC consistent with the triggered bunch; 

• energy sharing among the towers was used to reject electromagnetic showers from 
high-energy photons. 

These cuts, similar to those used in the previous ZEUS measurements [2,5], selected long- 
lived neutral hadrons which had not interacted with material before reaching the FNC. 
The sample was predominantly neutrons, with a small component of hadrons. The 
MC models Rapgap with pion exchange [24] and Lepto with soft color interactions [25] 
predict that the contribution is less than 2% above xl = 0.6, and increases slowly to 
approximately 20% at xl = 0.2. 

The above selection was used for the x^ measurements. For results requiring also a 
measurement of p^, the following additional requirements were imposed to ensure a well 
reconstructed position measurement in the FNT: 

• the highest pulse-height channel in each of the hodoscope planes was above the pedestal 
level, to select neutrons which showered before the FNT plane; 
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• shower profiles with more than one peak were rejected, to eliminate mismeasurement 
from shower fluctuations. 

The fraction of clean FNC neutrons passing the FNT cuts determined the FNT efficiency 
as a function of xl- The efficiency rises with neutron energy from 35% at xl = 0.2 to 
63% at Xl = 0.85, corresponding to the fraction of neutrons that shower before the FNT. 



3.2 Normalization 

The cross sections for leading neutron production presented here, (Tln? were normalized 
to the inclusive cross sections without a leading- neutron requirement, (Tine, as: 

Variations of this relative neutron yield, tln, with W"^ or indicate differences in the 
neutron-production mechanism. The acceptance for detecting different types of events in 
the central ZEUS detector in a small kinematic region of {W'^,Q'^) is independent of the 
neutron requirement; the acceptance cancels in the yield tln, so that: 

Here A^inc is the number of inclusive events in the sample and A^ln is the number of these 
events with a neutron tag, corrected for the acceptance of the forward neutron detectors. 
The acceptance of the central ZEUS detector varies with {W'^, Q'^)', if the neutron-tagged 
and inclusive events have different kinematic dependences, their acceptances integrated 
over a given {W"^, Q"^) region would be different. The mild violation of vertex factorization 
observed in the previous ZEUS measurement [2] indicates that such differences in accep- 
tance are less than 2%, and so were ignored. Thus the acceptance of the central ZEUS 
detector and associated systematic uncertainties do not affect the neutron yield. Only the 
acceptance of the forward neutron detectors together with its systematic uncertainties are 
relevant for measuring tln- 

For the DIS sample, a set of inclusive events was collected simultaneously and used to 
normalize the neutron data. For the inclusive photoproduction events, 49nb~^ of data 
were collected in a special run in 1996 when the proton beam energy was 820 GeV, for 
measurement of the photon-proton total cross section [23]. The normalization of the 2000 
photoproduction data was determined by assuming that the ratio of the photoproduction 
and DIS neutron yields, for any given neutron {x^^p^) kinematic region, was the same at 
the two proton beam energies: 



LN 



^DIS 

'ln 



^7P 
LN 



».DIS 

920 GeV ' LN 
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In this equation rL^(820GeV) was measured in the 1996 inclusive photoproduction sam- 
ple, rLN^(820 GeV) was measured in inclusive DIS data from the same running period and 
rL^^(920 GeV) was measured in the inclusive DIS data from the 2000 running period. The 
neutron kinematic region for the measured yields was xl > 0.2 and p"^ < 0.476 x| GeV^. 
This normalization procedure resulted in an uncertainty on the neutron yield in photopro- 
duction of 5.1%, predominantly from the limited statistics of the 1996 photoproduction 
data. 

3.3 Beamline and forward-detector simulation 

The acceptance of the forward neutron detectors, for the measurement of neutron yields, 
was determined from a simple one-particle MC simulation. The simulation accounts for 
the aperture and dead material along the neutron flight path, the measured proton beam 
position and pr spread, and the measured detector resolutions. 

Figure [5] shows a scatter plot of reconstructed hits in the FNT from a sample of DIS 
events. The irregular curve is the aperture expected from the MC simulation. Numerous 
events are reconstructed outside of the aperture, as the aperture is not a sharp boundary 
as modeled, but presents a varying amount of dead material over several millimeters 
transverse to the neutron flight path. The effect of this on the measured neutron yield is 
less than 2% and was ignored. 

The simulation also modeled significant amounts of dead material along the neutron flight 
path, primarily from stations S5 and S6 of the ZEUS leading proton spectrometer (LPS) 
shown in Fig. [2l The LPS was a set of Roman pots used to measure protons scattered 
at very low angles [26]. The elements of these stations were measured after the LPS 
was removed from the HERA tunnel in 2000 and implemented in the simulation. The 
positions of these elements in the simulation were adjusted to reproduce the data. For 
example, the deficit of events observed near Y ~ 10.5 cm, seen clearly in the vertical slice 
of the scatter plot in Fig. [6l determined the vertical alignment of the LPS. There is good 
agreement between the simulation and the data distributions. 

During operation, the LPS was normally in one of two positions: extracted, or inserted for 
data taking. Separate dead material maps were made for the two positions. Data collected 
during periods when the LPS was moving were rejected. The results were determined 
separately for the two LPS positions and combined according to the luminosity taken in 
each position. The difference with the results obtained for each position was taken as a 
measure of the systematic uncertainty from the dead material map. 

The reconstruction of the neutron scattering angle, 6n, requires knowledge of the zero- 
degree point. This was determined by generating a symmetric distribution of neutrons. 
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passing it through the simulation, and fitting the reconstructed distribution to the data. 
An example of such a fit for the zero-degree position in the vertical plane is shown in Fig. [61 
Considering different input distributions, and taking into account the uncertainties in the 
dead material map, the beam zero-degree point was determined to an accuracy of ±0.2 cm 
in both X and Y . 

The simulation also takes into account the energy resolution of the FNC, the position 
resolution of the FNT, and the inherent p^-spread of the HERA proton beam. The 
latter was measured in the reaction 7p pp, with the p decay products measured in the 
ZEUS central detector and the final-state proton measured with the LPS [26]. The beam 
PT-spread corresponds to a smearing of the zero-degree point by 0.45 cm horizontally 
and 1.0 cm vertically, significantly larger than the FNT resolution. The spectrum of 
generated neutrons was tuned to match the xl and p"^ distributions separately for the 
DIS and photoproduction samples. 

These MC distributions were then used to correct the data for all acceptance and smearing 
effects. For distributions requiring a position measurement, the correction for the x^- 
dependent FNT efficiency was also applied. 

3.4 Systematic uncertainties 

The dominant effects contributing to the systematic uncertainties arose from: 

• the beam zero-degree point; 

• the dead material map; 

• the FNC energy scale; 

• the distributions, input for ^/.-distribution acceptance. 

The systematic uncertainties were typically 5-10% of the measured quantities, and are 
shown as shaded bands in the figures. The variation of the energy scale shifted the points 
in xl- The other systematic variations amount to a change in acceptance resulting in a 
correlated shift of the neutron yields. The small uncertainties from the assumption of 
acceptance cancellation in tln and from the aperture edge were ignored. 

Corrections for efficiency of the cuts and backgrounds were applied to the normalization 
of the neutron yields. The corrections, similar to those of the previous ZEUS measure- 
ment [2], were: 

• false veto counter firing: +10.6 ± 1.0%, determined from randomly triggered events; 

• veto counter inefficiency: —2.5 ± 1.0%, determined from 2-mip distributions in the 
veto counter; 
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• backsplash from neutron showers: +1.5 ± 1.5%, determined from timing information 
in the veto counter giving the fraction of late time vetoes; 

• neutrons from proton beam-gas interactions: —1.4 ±0.3%, determined from randomly 
triggered events. 

The overall systematic uncertainty on the normalization of the neutron yield from these 
corrections, not included in the shaded bands of the figures, was ±2.1%. The yield in pho- 
toproduction has an additional uncertainty of ±5.1% from the normalization procedure 
described in Section 13. 2[ 

4 Results 

4.1 Neutron xl and distributions in DIS 

Figure [7] shows the normalized differential distribution {l/ai^c)dcri^^/dxL for neutrons 
in DIS with scattering angles 6n < 0.75 mrad, corresponding to the kinematic range 
pI < 0.476 x|GeVl It rises from the lowest xl due to the increase in p'^ space, reaches 
a maximum near xl = 0.7, and falls to zero at the endpoint xl = 1. These results are 
consistent with the previous ZEUS measurement [2]. Integrating this distribution, the 
total leading-neutron yield for the measured region is: 



rMQ"^ > 2GeV^XL > 0.2, < 0.476 x^GeV^) = 0.0885 ± 0.0002 (stat.) 1°;°°^^ (sys.). 



Here the systematic uncertainty includes the overall ±2.1% scale uncertainty. 

The corrected p"^ distributions in DIS for different xl bins are shown in Fig. [8] and sum- 
marized in Table [TJ They are presented as normalized doubly differential distributions 
{l/ainc)d'^o'i,K/dxLdp'^. The bins in p'^ are at least as large as the resolution, which is 
dominated by the pt spread of the proton beam. The varying p"^ ranges of the plots 
are due to the aperture limitation. The line on each plot is a fit to the functional form 
daui/dpip oc exp(— The distributions are compatible with a single exponential within 
the statistical and uncorrelated systematic uncertainties. Thus, with the parameterization 



the neutron {xl.Pt) distribution is fully characterized by the slopes h^xi) and intercepts 
o-^xl) = (l/crinc)'i^o"LN/'^a^L'^PTlp2,=o- '^^^ dependences of the intercepts and the slopes on 
Xl are shown in Figs. [9] and [lO] and summarized in Table [21 The systematic uncertainties 
were evaluated by making the variations listed in Section 13.41 and repeating the fits. The 
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intercepts fall rapidly from the lowest x^, are roughly constant in the region xl = 0.45- 
0.75, and fall to zero at the endpoint xl = 1. Below xl = 0.32, the slopes are consistent 
with zero and are not plotted; they rise linearly in the range 0.3 < Xl < 0.85 to a value 
of 6 ~ 8 GeV~^, and then decrease slightly at higher xl- Figure [TT] shows the ratio a/h for 
the region 0.32 < xl < 1 where 6 > 0, which can be taken as the leading-neutron yield 
integrated over values from zero to infinity, assuming that the p\ distributions remain 
an exponential also beyond the measured range. This distribution, integrated over xl 
in the range 0.32 < < 1, corresponds to a yield: 

rLN(<3^ > 2GeV^XL > 0.32) = 0.159 ±0.008 (stat.) (sys.). 

Here the systematic uncertainty includes the overall 2.1% scale uncertainty. 

4.2 dependence of leading neutron production 

To investigate the dependence of leading-neutron production, the full DIS sample 
was divided into three subsamples depending on the Qda range, with the additional cut 
yjB > 0.02. The kinematic regions in and W for all DIS and photoproduction samples 
are summarized in Table [31 

Figure [I2] shows the xl distributions and Fig. [13] the pf. slopes in the range p\ < 
0.476 x\ GeV^ for the photoproduction and three DIS subsamples. There is a trend of in- 
creasing neutron yield with increasing Q^, a clear violation of vertex factorization. There 
is a large increase between the photoproduction region and the low- and mid-Q^ DIS 
regions, in which the data are nearly independent. There is then a smaller, but sig- 
nificant, increase between the mid-Q^ and high-Q^ regions. The effect of the different W 
ranges for the DIS and photoproduction samples on the neutron yield was less than 5%, 
as evaluated by restricting the DIS sample to low- and high-iy regions. The slopes for 
all three DIS samples are equal within the uncertainties. The slopes for photoproduction 
are higher in the region 0.6 < xl < 0.9. 

The total neutron yields integrated over xl > 0.2 for the four samples are summarized in 
Table[ll The doubly differential distributions 1/ ai^icdcTLN / dxidp^ for the photoproduction 
and three DIS subsamples are summarized in Table [5l and the intercepts and slopes of 
the exponential fits are summarized in Table [61 

To investigate the differences between the photoproduction and DIS regimes further, the 
effects of energy calibration and beam position drifts were minimized by using only the 
subset of DIS data collected simultaneously with the photoproduction data. The DIS 
sample without ?/jb or Qda ^^^^ used to maximize the statistical precision of the 
comparison. 
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The ratio of the normahzed differential distributions 




1 da^^ 




for the region p\ < 0.476 x|^GeV^ is shown in Fig. (TH In the range 0.2 < xl < 0.4, the 
ratio drops shghtly but rises for higher xl values, exceeding unity for xl > 0.9. The 
deviation of the ratio from unity is a clear violation of vertex factorization. The ratio of 
the intercepts for photoproduction and DIS, which has a nearly identical behavior to that 
of Fig. [131 is not shown. 

The distributions for both samples, normalized to unity at p"^ = 0, are shown in Fig. [151 
The photoproduction distributions are steeper in the range 0.6 < xl < 0.9, with relatively 
fewer neutrons at high The difference of the slopes Ab = h{Q'^ < 0.02 GeV^) — h^Q"^ > 
2 GeV^) is less sensitive to systematic effects than each of the individual slopes. These 
values are plotted in Fig. [161 The slopes for photoproduction are larger in the range 
0.6 < Xl < 0.9, with Ab = 0.5-1.0 GeV~^, qualitatively consistent with the violation of 
vertex factorization expected from absorption as discussed in the introduction. 

4.3 Comparison to leading protons 

ZEUS has also reported xl distributions for leading protons in the kinematic region 
< 0.04 GeV^ [4]. The neutron xl distributions were also measured in the same region, 
using the FNT measurement of p"^. The results for DIS are compared in Fig. [T71 There 
are approximately twice as many protons as neutrons in the range 0.6<xl<0.9. If 
only isospin-1 particle exchanges contributed to proton production, there should be half 
as many protons as neutrons. Thus, exchanges of particles with different isospins such as 
isoscalars must be invoked to account for the observed proton rate [11]. 

The slopes of the p"^ distributions for leading protons and neutrons in DIS are shown in 
Fig. [m Note that the p^ range for the proton measurement, p"^ < 0.5 GeV^, is larger than 
for the neutron measurement. The two samples have similar values of b near xl ~ 0.8, 
the region where pion exchange is expected to dominate for both processes [11]. 

5 Comparison to models 

In this section the data are compared to several models. First the data are compared 
to various MC models for the simulation of DIS events. A comparison is then made 
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to models incorporating only pion exchange. Next, more sophisticated models including 
the effects of absorption of the neutron are considered. Finally, a model incorporating 
enhanced absorption with pion and additional secondary Regge exchanges is discussed. 

5.1 Monte Carlo models 

Most MC models generate leading neutrons from the fragmentation of the proton rem- 
nant [27]. Some models also incorporate additional processes to simulate diffraction and 
leading baryon production. The leading neutron x l distribution, intercepts and slopes in 
DIS are compared to two MC models in Fig. [191 The models presented here are Rap- 
gap [24] and Lepto [28]. The proton PDF parameterization used was CTEQ5L [29]. 
With only standard proton-remnant fragmentation, the models are lower than the data 
in the normalization of the xl distribution and intercepts, and are peaked at lower xl- 
They do not show the observed xl dependence of the slopes. Other models incorporating 
only standard fragmentation, Ariadne [30] and Cascade [31] for DIS, and Pythia [32] 
and Phojet [33] for photoproduction, give a similarly poor description of the data. 

Lepto has the option to implement soft color interactions (SCI) [25] to produce rapidity 
gaps observed in diffractive events. This model gives a qualitative description of the 
leading proton xl distribution [4], including the diffractive peak, although it predicts 
too few protons in the central xl region. The predictions for leading neutrons, with the 
probability of SCI set to 0.5, are shown in Fig. ^\ Lepto SCI comes close to the data in 
the shape and normalization of the xl distribution. The intercepts also exhibit a shoulder 
in the distribution near xl ~ 0.8 similar to that in the data. It does not, however, give 
the observed strong xl dependence of the slopes. 

Rapgap also includes Pomeron exchange to simulate diffractive events, and pion ex- 
change to simulate leading baryon production. These processes are mixed with standard 
fragmentation according to their respective cross sections. The PDF parameterizations 
used here were the HI LO fit 2 [34] for the Pomeron and GRV-P LO fit [35] for the 
pion. The predictions for leading neutrons from a mixture of these exchanges and stan- 
dard fragmentation are also shown in Fig. [191 The model well reproduces the shape of 
the Xl distribution and intercepts, although it predicts more neutrons than are observed. 
The model also shows the strong xl dependence of the slopes in the data, although the 
predicted values of the slopes are systematically larger than the data. 

5.2 Pure pion exchange 

In the Regge factorization relation discussed in the introduction, the leading neutron xl 
distribution is a product of the pion flux, fn/p, and the 7*7r cross section, a^*-,^. However, 
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if (T^*^ is assumed to be independent of t, the pt distribution of the produced neutrons is 
determined only by the pion flux /tt/p- The slopes can be compared to various parame- 
terizations of the flux. Although /^r/p is not an exponential in p"^, at fixed xl the models 
can be fit to the form exp{—bp^) using the same binning as the data, and the resulting 
K^l) values compared to the measurements. All of the parameterizations in the litera- 
ture [6-9, 11] give values for the slopes larger than the data. Most of them also have the 
wrong xl dependence of the slopes. The models that most resemble the data are shown 
in Fig. [201 The simple model of Bishari [6], with the form factor F{xL^t) = 1, is closest 
in magnitude to the data. Other models with more detailed parameterizations show the 
turnover of the slopes near xl ~ 0.85 [7]. The model of Holtmann et al. [8] is used for 
pion exchange in the Rapgap-tt MC. The values for b from Rapgap-tt in Fig. [12] are 
lower than this curve because Rapgap also includes a contribution from standard frag- 
mentation, which has flatter p"^ distributions than pion exchange. None of these models, 
based on pion exchange alone, describes the data. 

5.3 Pion exchange with neutron absorption 

As discussed in the introduction, in a geometrical picture, neutron absorption may occur 
for large photon sizes and small n-vr separations. The former is inversely related to Q^, and 
so is largest for photoproduction; thus more absorption is expected in photoproduction 
than in DIS. The n-ir separation is the Fourier conjugate oi pt, and the distribution 
of Tmr is glveu by the Fourier transform of fn/piPT)- Parameterizations of the pion flux in 
general show that the mean value of increases with x/,, so more absorption is expected 
at lower xl than at higher xl- The dashed curve in Fig. [21] is the expectation for the 
suppression of leading neutrons in photoproduction relative to DIS from a model of pion 
exchange with neutron absorption [14]. Although the curve lies below the data, it follows 
the same trend. The 7p interaction has a power- law dependence a oc W"^^, with different 
values of A for DIS and photoproduction. Assuming that 77r interactions have the same 
dependence, and recalling that W^^r = ^/l — xiW^p, the ratio of photoproduction and 
DIS cross sections is proportional to (1 — xl)^^. Previous ZEUS measurements of A in 
photoproduction [23] and DIS [36,37] give AA ~ —0.13. Applying this to the absorption 
suppression factor results in the solid curve in Fig. [211 Within the normalization uncer- 
tainty of 5.1%, the data are well described by the absorption model with this correction 
for different W dependences. Hence such a geometric absorption model can account for 
the differences between the xl distributions in DIS and photoproduction. 

Also shown in Fig. [21] is the NSZ model [13] which employs the optical theorem to- 
gether with multi-Pomeron exchanges to describe all possible rescattering processes of 
the leading hadron, resulting in absorptive effects. With the correction for different W 
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dependences, the prediction is close in magnitude to the data, but does not have as steep 
an xl dependence. 

5.4 Enhanced neutron absorption and secondary exchanges 

Recently a new calculation (KKMR) of pion exchange with neutron absorption based 
on multi-Pomeron exchanges has become available [15]. The pion exchange is based on 
the Bishari flux. In addition to the rescatterings implemented in the earlier model [13], a 
small contribution from rescattering on intermediate partons in the central rapidity region 
is also included. The model also accounts for the migration of neutrons in {xi^p^) after 
rescattering. The prediction for the neutron xl distribution for photoproduction, where 
rescattering is most important, is shown by the dashed curve in Fig. [22j The model gives 
a fair description of both the shape and normalization of the data. The loss of neutrons 
through absorption is approximately 50%; this is consistent with the deviation from the 
cr^^/cr^p = 2/3 prediction of the additive quark model that was noted in the previous 
ZEUS measurement [2]. Within this model, the present data can be used to constrain the 
gap-survival probability, one of the crucial inputs to calculations of diffractive interactions 
at the LHC — both hard, such as central exclusive Higgs production, and soft, such as 
those giving rise to the diffractive pile-up events [38] . The prediction of this model for the 
slopes h in DIS is shown in Fig. [231 As for the pure pion-exchange calculations, the model 
predicts larger values of h than seen in the data. This model does give a fair prediction 
for the magnitude of the difference of the slopes in photoproduction and DIS, as shown 
in Fig. [221 

More recently this model was extended to include, in addition to pion exchange, the 
exchange of secondary (p, 02) Reggeons [16]. This extended model gives a fair description 
of the shape and normalization of the xl distribution in photoproduction, as seen in 
Fig. [22I Since there are additional exchanges the model gives a prediction for the xl 
distribution which is higher than for pion-exchange alone. As shown in Fig. [231 the model 
with secondary exchanges also gives a good prediction for the slopes. Its description of 
the slope differences is also close in magnitude to the data, as seen in Fig. [231 

6 Summary 

The Xl and p\. distributions of leading neutrons in photoproduction and DIS events at 
HERA have been measured. The x^ distributions for the measured region On < 0.75 mrad 
rise from the lowest xl due to the increasing pt phase space, reach a maximum near 
Xl = 0.7, and fall to zero at the endpoint = 1. The p'^ distributions are well described 
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by an exponential da/dp^ = aexp(— The intercepts a fall rapidly from the lowest 
xl, are roughly constant in the region xl = 0.45-0.75, and fall to zero at the endpoint 
xl = 1. The exponential slopes b rise linearly with xl in the range 0.3 < xl < 0.85 to a 
value of 6 ~ 8 GeV~^, and then decrease slightly at higher xl- 

The neutron yield rises monotonically with from the photoproduction region < 
0.02 GeV^ to the high-Q^ dIS region g2 > 20GeV^ The relative rise in yield is greatest 
near xl = 0.5 and becomes less significant at higher xl- The slopes of the distribu- 
tions do not change significantly within the DIS region > 2 GeV^, but the slopes in 
photoproduction exhibit a small increase over those in DIS for 0.6 < xl < 0.9. 

In the kinematic region 0.6 < xl < 0.9 and < 0.04 GeV^ there are approximately 
twice as many leading protons as neutrons. This indicates that leading proton production 
proceeds through exchanges in addition to pure isovector (e.g. pion) exchange. The slopes 
b for leading protons agree with the neutron slopes near xl ~ 0.8 where pion exchange is 
expected to dominate both processes. 

Monte Carlo models commonly in use for the simulation of DIS and photoproduction 
events which implement standard fragmentation of the proton remnant do not describe 
the leading-neutron data. They predict fewer neutrons, concentrated at lower x^. They 
also predict smaller p"^ slopes and do not have the strong xl dependence of the data. The 
inclusion of soft color interactions gives a reasonable description of the xl distributions, 
but again fails to predict the p'^ slopes. A mixture of processes including standard frag- 
mentation, diffraction and pion exchange gives a good description of the xl distributions 
and the xl dependence of the slopes, although they are larger than the data. 

The measured b{xL) dependence in DIS has been compared with various pion-exchange 
models. All models give values larger than the data. The simplest model is closest in 
magnitude to the data; other models reproduce the measured shape of b^xi). 

The dependence of the neutron yield and p^ slopes is consistent with absorption 
models where neutrons from pion exchange with smaller n-vr separations are lost through 
rescattering on larger photons. The photon size increases with decreasing Q^, and the 
mean n-vr separation is smaller in the mid-x^ range than at higher xl- The result is 
a depletion of neutrons with decreasing Q^, with the depletion greater at mid-x^ than 
at higher xl-, as seen in the data. The loss of neutrons with small n-ir separations, 
corresponding to large pf,, also explains the larger pip slopes measured in photoproduction 
than in DIS. 

A model of neutron production through pion exchange, incorporating enhanced neutron 
absorption and migration of the neutrons in {xl^Pt) after rescattering, gives a fair descrip- 
tion of the shape and normalization of the x^ distributions in DIS and photoproduction, 
and of the difference in the p\ slopes b between the two sets. However, as with pure 
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pion exchange, it predicts too high a value for b. Extending the model to include also p 
and 02 exchanges still gives a fair description of the shape and normalization of the xl 
distributions in DIS and photoproduction, and also good descriptions of the p"^ slopes and 
of the differences between the two sets. 
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(xt) 




1/(T;„^ d ai K\ I dx T driin ( GeV 


0.20-0.32 


0.27 


7.96 -10-^ 
2.51 -10-3 
4.86 -lO-^ 
7.95 -10-3 
1.19 -10-2 
1.65 ■10~2 


2.031 ± 0.064 
1.881 ± 0.064 
1.955 ± 0.070 
2.002 ± 0.068 
1.962 ± 0.063 
2 001 ± 062 


0.32-0.42 


0.37 


1.99 -10-3 
6.46 -10-3 
1.24 -10-2 
2.04 -10-2 
3.03 -10-2 
4.23 ■10~2 


1.541 ± 0.035 
1.454 ± 0.038 
1.445 ± 0.039 
1.527 ± 0.037 
1.440 ± 0.034 
1 256 ± 031 


0.42-0.50 


0.46 


3.42 -10-3 
1.11 ■10~2 
2.14 -10-2 
3.51 -10-2 
5.22 -10-2 
7.28 ■10"2 


1.336 ± 0.026 
1.211 ± 0.029 
1.218 ± 0.029 
1.217 ± 0.027 
1.147 ± 0.025 
1 005 ± 023 


0.50-0.54 


0.52 


4.84 -10-2 
1.58 -10-2 
3.03 -10-2 
4.97 ■10~2 
7.40 -10-2 
1.03 -10"^ 


1.274 ± 0.030 
1.218 ± 0.035 
1.200 ± 0.034 
1.110 ± 0.031 
0.946 ± 0.026 
836 ± 026 


0.54-0.58 


0.56 


5.64 -10-3 
1.84 -10-2 
3.53 ■10"2 
5.80 -10-2 
8.62 ■10"2 
1.20 -10-^ 


1.328 ± 0.028 
1.144 ± 0.031 
1.166 ± 0.031 
1.016 ± 0.027 
0.904 ± 0.024 
0.767 ± 0.023 



Table 1: The normalized doubly differential distributions {l/ai^c)d'^o'L^i/dxLdpj 
for the full DIS sample. Only statistical uncertainties are shown. 
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Table [U (cont.) 



xl range 




P^GeV^) 


l/cTinc daiM/dxLdpT ( GeV~^) 


0.58-0.62 


0.60 


6.50 -10-3 
2.12 -10-2 
4.07 -10-2 
6.68 -10-2 
9.94 -10"^ 
1.39 -10-1 


1.286 ± 0.025 
1.175 ± 0.029 
1.112 ± 0.028 
0.976 ± 0.024 
795 ± 020 
0.654 ± 0.019 


0.62-0.66 


0.64 


7.42 -10-3 
2.42 -10-2 
4.65 -10-2 
7.63 -10-2 
1.14 -10"^ 
1.58 -10-1 


1.296 ± 0.024 
1.157 ± 0.027 
1.034 ± 0.025 
0.926 ± 0.022 
703 ± 018 
0.571 ± 0.017 


0.66-0.70 


0.68 


8.39 -10-3 
2.74 -10-2 
5.27 -10^2 
8.64 -10-2 
1.29 -10"^ 
1.79 -10-1 


1.297 ± 0.022 
1.112 ± 0.025 
0.990 ± 0.022 
0.821 ± 0.019 
581 ± 015 
0.456 ± 0.014 


0.70-0.74 


0.72 


9.42 -10-3 
3.08 -10-2 
5.92 -10-2 
9.71 -10-2 
1.44 -10"^ 
2.01 -10-1 


1.306 ± 0.021 
1.126 ± 0.023 
0.914 ± 0.020 
0.707 ± 0.016 
0.524 ± 0.013 
0.388 ± 0.012 


0.74-0.78 


0.76 


1.05 -10-2 
3.44 -10-2 
6.61 -10-2 
1.08 -10-1 
1.61 -10-1 
2.25 -10-1 


1.280 ± 0.019 
1.030 ± 0.021 
0.837 ± 0.018 
0.617 ± 0.014 
0.403 ± 0.011 
0.308 ± 0.010 
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Table [I](cont.) 



Xl range 






l/cTinc daui/dxidpT ( GeV"^) 


0.78-0.82 


0.80 


1.16 -10-2 
3.81 -10-2 
7.33 -lO-^ 
1.20 -10-1 
1.79 -10"^ 
2.49 -10-^ 


1.180 ± 0.018 
0.920 ± 0.019 
0.714 ± 0.016 
0.482 ± 0.012 
304 ± 009 
0.213 ± 0.008 


0.82-0.86 


0.84 


1.28 -10-2 
4.21 -lO-^ 
8.10 -10-2 
1.33 -10"^ 
1.98 -10"^ 
2.75 -10-1 


1.000 ± 0.016 
0.783 ± 0.017 
0.568 ± 0.014 
0.374 ± 0.010 
208 ± 007 
0.139 ± 0.006 


0.86-0.90 


0.88 


1.41 -10-2 
4.63 -10-2 
8.90 -10-2 
1.46 -10-^ 
2.17 -10"^ 
3.03 -10-^ 


0.719 ± 0.013 
0.537 ± 0.013 
0.363 ± 0.010 
0.232 ± 0.007 
135 ± 005 
0.082 ± 0.004 


0.90-0.95 


0.92 


1.54 -10-2 
5.07 -10-2 
9.74 -10-2 
1.60 -10-^ 
2.38 -10"^ 
3.31 -10-^ 


0.347 ± 0.007 
0.267 ± 0.008 
0.172 ± 0.006 
0.110 ± 0.004 
0.059 ± 0.003 
0.034 ± 0.002 


0.95-1.00 


0.97 


1.72 -10-2 

5.64 -10-2 
1.08 -10-^ 
1.78 -lO^i 

2.65 -10-^ 
3.69 -10-^ 


0.054 ± 0.002 
0.043 ± 0.002 
0.032 ± 0.002 
0.015 ± 0.001 
0.009 ± 0.001 
0.005 ± 0.001 
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xl range 






a 


( GeV-2) 


h ( GeV-2) 





20- 


-0 


32 





27 


1 


958 


±0 


045l°:J^ 


-0.94 ±2.401?:}^ 





32- 


-0 


42 





37 


1 


551 


±0 


not; +0.10 
UZO_Q08 


3.93 ±0 


co+om 

O°-0.90 





42- 


-0 


50 





46 


1 


325 


±0 


Q2n+0.10 


3.38 ±0 


07+0. 67 
' -0.54 





50- 


-0 


54 





52 


1 


324 


±0 


Q24+0.12 


4.32 ± 


oq+0.68 
'^'5-0.53 





54- 


-0 


58 





56 


1 


327 


±0 


Q22+0.10 


4.53 ± 


97+0.45 
^ ' -0.28 





58- 


-0 


62 





60 


1 


336 


±0 


Q21+0.14 

u^J--o.io 


5.09 ±0 


99+0.55 
^^-0.60 





62- 


-0 


66 





64 


1 


343 


±0 


Q2n+o.i6 

U^U_o.o8 


5.44 ± 


1 q+0.82 





66- 


-0 


70 





68 


1 


357 


±0 




6.23 ±0 


1 7+0.50 
' -0.45 





70- 


-0 


74 





72 


1 


370 


±0 


0181°:!^ 


6.53 ±0 


1 c,+0.72 
-■-^-0.58 





74- 


-0 


78 





76 


1 


344 


±0 


0171°:°? 


7.06 ± 


1/1+0.59 
-■-^-0.61 





78- 


-0 


82 





80 


1 


258 


±0 


UJ- ' -0.04 


7.66 ± 


14+0.73 
-■-^-0.57 





82- 


-0 


86 





84 


1 


098 


±0 


QIC- +0.05 
'-'-'-^-0.07 


8.04 ±0 


1/1+0.65 
-'-'^-0.54 





86- 


-0 


90 





88 





782 


±0 


niQ+o.09 


8.03 ±0 


1 c:+0.46 
-■-^-0.45 





90- 


-0 


95 





92 





387 


±0 




7.79 ± 


1 7+0.60 
' -0.45 





95- 


-1 


00 





97 





063 


±0 


0021°:°! 


7.46 ± 


90+0.59 
'^°-0.47 



Table 2: The intercepts a and slopes b from the exponential parameterization 
of the differential cross defined in Section for the full DIS sample. Statistical 
uncertainties are listed first, followed by systematic uncertainties, not including 
an overall normalization uncertainty of 2.1% on the intercepts. The systematic 
uncertainties are largely correlated between all points. 
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sample 


g2 range ( GeV^) 


(g2) ( GeV^) 


ly range ( GeV) 


{W) ( GeV) 


IP 


g2 < 0.02 


4 X 10-^ 


150 < < 270 


215 


full DIS 


g' >2 


13 


W < 250 


95 


low-g2 DIS 


2 < gi)A < 5 


2.7 


50 < < 250 


140 


mid-g2 DIS 


5 < QIj^ < 20 


8.9 


50 < < 250 


132 


high-g2 DIS 


20 < < 120 


40 


50 < 1^ < 250 


131 



Table 3: Kinematic regions for each of the data samples. The 7p ranges were 
estimated from a simulation of photoproduction including the positron tagger used 
to trigger these data. The ranges for the low-, mid-, and high-Q^ DIS samples 
are the limits on the double angle variable Q'jjj^. All other ranges and means for DIS 
were estimated using the DJANGOH 1.1 [39] generator, where the CTEQ4D [40] 
parton- density parameterizations were used. 



sample 




IP 


0.0700 ±0.0004 +°:°^ 


full DIS 


0.0885 ±0.0002 +°:°°^9 


low-g2 DIS 


0.0837 ±0.0003 t^^^ 


mid-g2 DIS 


0.0843 ±0.0003 t^^^^ 


high-g2 DIS 


0.0913 ±0.0006 t^^^ 



Table 4: tln in the region xl > 0.2, < 0.476x| GeV'^ , for the photoproduction 
and DIS samples. Statistical uncertainties are listed first, followed by systematic 
uncertainties, not including an overall normalization uncertainty of 2. 1 %. 



24 



X]^ range 






l/(Tinc rfcTLN / dxLdpx ( GeV 2) 


IP 


10W-Q2 DIS 


inid-Q2 DIS 


high-Q2 DIS 


0.20-0.50 


0.38 


8.40 •10"'' 
2.43 -10-3 
4.86 -10-3 
7.94 -10-3 
1.19 •lO"^ 
1.65 -10-2 


1.501 ± 0.108 
1.289 ± 0.090 
1.127 ± 0.076 
1.232 ± 0.084 
1.064 ± 0.075 
1.021 ± 0.068 


1.491 ± 0.058 
1.342 ± 0.052 
1.370 ± 0.052 
1.417 ± 0.054 
1.359 ± 0.053 
1.374 ± 0.052 


1.721 ± 0.063 
1.500 ± 0.056 
1.451 ± 0.053 
1.528 ± 0.057 
1.451 ± 0.055 
1.589 ± 0.056 


1.848 ± 0.115 
1.749 ± 0.106 
1.757 ± 0.104 
1.809 ± 0.110 
1.803 ± 0.109 
1.805 ± 0.106 


0.50-0.58 


0.54 


4.84 -10-3 
1.58 -10-2 
3.03 -10-2 
4.97 •lO"^ 
7.39 •lO"^ 

1.03 -10-^ 


0.988 ± 0.056 
0.791 ± 0.055 
0.812 ± 0.058 
0.784 ± 0.052 
0.654 ± 0.043 
0.591 ± 0.041 


1.196 ± 0.036 
1.077 ± 0.042 
1.027 ± 0.040 
1.005 ± 0.037 
0.880 ± 0.033 
0.768 ± 0.031 


1.247 ± 0.037 
1.096 ± 0.042 
1.042 ± 0.041 
1.009 ± 0.037 
0.857 ± 0.032 
0.760 ± 0.031 


1.452 ± 0.071 
1.220 ± 0.079 
1.191 ± 0.078 
1.187 ± 0.072 
0.934 ± 0.060 
0.782 ± 0.056 


0.58-0.66 


0.62 


6.50 -10-3 
2.12 -10-2 
4.08 -10-2 
6.68 -lO-^ 
9.94 •lO"^ 
1.39 -10-1 


0.991 ± 0.047 
0.849 ± 0.051 
0.892 ± 0.055 
0.753 ± 0.045 
0.563 ± 0.035 
0.451 ± 0.031 


1.192 ± 0.031 
1.062 ± 0.035 
1.025 ± 0.034 
0.906 ± 0.030 
0.714 ± 0.025 
0.575 ± 0.023 


1.261 ± 0.032 
1.122 ± 0.037 
1.068 ± 0.035 
0.889 ± 0.030 
0.748 ± 0.025 
0.605 ± 0.023 


1.221 ± 0.055 
1.131 ± 0.065 
1.151 ± 0.065 
1.090 ± 0.059 
0.837 ± 0.048 
0.626 ± 0.042 


0.66-0.74 


0.70 


8.38 -10-3 
2.74 -10-2 

5.27 -10-2 
8.64 -10-2 

1.28 -10-^ 
1.79 -10-1 


1.083 ± 0.044 
0.882 ± 0.047 
0.815 ± 0.045 
0.688 ± 0.039 
0.451 ± 0.027 
0.305 ± 0.021 


1.303 ± 0.028 
1.116 ± 0.032 
0.935 ± 0.028 
0.751 ± 0.024 
0.530 ± 0.018 
0.426 ± 0.017 


1.263 ± 0.028 
1.093 ± 0.031 
0.961 ± 0.029 
0.739 ± 0.023 
0.560 ± 0.019 
0.428 ± 0.017 


1.309 ± 0.050 
1.104 ± 0.056 
1.035 ± 0.054 
0.745 ± 0.042 
0.532 ± 0.033 
0.448 ± 0.031 



Table 5: The normalized doubly differential distributions {l/a'inc)d'^cri^N/dxLdp'^ 
for the photoproduction and three DIS samples. Only statistical uncertainties are 
shown. 



25 



Table [5](cont.) 



XL range 




p|(GeV2) 


l/<7inc dau^/dxLdp^ ( GcV ^) 


7P 


1ow-Q2 dIS 


mid-Q2 DIS 


high-Q2 DIS 


0.74-0.82 


0.78 


1.05 -10-2 
3.44 -10-2 
6.61 -10-2 
1.08 -IQ-^ 
1.61 -10-1 
2.25 -10-1 


1.111 ± 0.041 
0.884 ± 0.043 
0.769 ± 0.042 
0.494 ± 0.028 
0.285 ± 0.018 
0.209 ± 0.016 


1.250 ± 0.025 
1.027 ± 0.027 
0.824 ± 0.024 
0.587 ± 0.018 
0.361 ± 0.013 
0.267 ± 0.012 


1.181 ± 0.024 
0.986 ± 0.027 
0.779 ± 0.023 
0.541 ± 0.018 
0.362 ± 0.013 
0.243 ± 0.011 


1.251 ± 0.044 
0.915 ± 0.046 
0.846 ± 0.043 
0.576 ± 0.033 
0.399 ± 0.025 
0.289 ± 0.022 


0.82-0.90 


0.86 


1.28 -10-2 
4.21 -10-2 
8.09 -10-2 
1.33 -lO^i 
1.97 -10-1 
2.75 -10-1 


0.876 ± 0.034 
0.688 ± 0.036 
0.538 ± 0.033 
0.293 ± 0.020 
0.134 ± 0.011 
0.098 ± 0.010 


0.899 ± 0.019 
0.732 ± 0.021 
0.528 ± 0.017 
0.335 ± 0.013 
0.195 ± 0.009 
0.116 ± 0.007 


0.894 ± 0.019 
0.671 ± 0.020 
0.477 ± 0.016 
0.329 ± 0.012 
0.170 ± 0.008 
0.113 ± 0.007 


0.838 ± 0.032 
0.647 ± 0.035 
0.477 ± 0.029 
0.268 ± 0.020 
0.174 ± 0.015 
0.116 ± 0.012 


0.90-1.00 


0.93 


1.54 -10-2 
5.06 -10-2 
9.74 -10-2 
1.60 -10-1 
2.38 -lO^i 
3.31 -lO"! 


0.242 ± 0.015 
0.166 ± 0.014 
0.130 ± 0.012 
0.072 ± 0.008 
0.039 ± 0.005 
0.022 ± 0.004 


0.236 ± 0.008 
0.180 ± 0.008 
0.117 ± 0.006 
0.077 ± 0.005 
0.041 ± 0.003 
0.019 ± 0.002 


0.228 ± 0.007 
0.156 ± 0.007 
0.115 ± 0.006 
0.065 ± 0.004 
0.035 ± 0.003 
0.023 ± 0.002 


0.201 ± 0.012 
0.142 ± 0.013 
0.093 ± 0.010 
0.057 ± 0.007 
0.034 ± 0.005 
0.019 ± 0.004 
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sample 


xl range 




a ( GeV"^) 


h ( GeV-2) 


IP 





20- 


-0 


50 





38 


1.386 ±0 


068t°:}° 


20.32 ± 5.271^:^^ 







50- 


-0 


58 





54 


0.949 ± 


04llHo^ 


4.72 ± 0.78l[j:^^ 







58- 


-0 


66 





62 


1.037 ±0 


0391°:^^ 


5.82 ± 0.531°:™ 







66- 


■0 


74 





70 


1.149 ±0 


038l°:ll 


7.13±0.40l^:?^ 







74- 


-0 


82 





78 


1.217±0 


0381°:^° 


8.32 ± 0.33l[j:l 







82- 


■0 


90 





86 


1.009 ±0 


'-''-'^-0.09 


9.23 ± QMtoM 







90- 


-1 


00 





93 


0.266 ±0 


0141H^ 


7.91 ± 0.44l°:^t 


iow-g2 





20- 


■0 


50 





38 


1.412 ±0 


038l°:J° 


2.21 ±2.92l|}^ 


DIS 





50- 


■0 


58 





54 


1.196 ±0 


0291°:°? 


4.20 ± 0.44l°:^J 







58- 


■0 


66 





62 


1.238 ±0 




5.40 ± QMtofo 







66- 


-0 


74 





70 


1.359 ±0 




6.88 ± 0.231°:^^ 







74- 


-0 


82 





78 


1.346 ±0 


r|9q+0.05 
"-''^'-'-0.06 


7.65 ± 0.19l°:^^ 







82- 


-0 


90 





86 


1.005 ±0 


0191°:°^ 


8.10±0.18l|]:^^ 







90- 


-1 


00 





93 


0.266 ±0 


oo8l°:J° 


7.94 ± 0.251°:^^ 


mid-Q^ 





20- 


-0 


50 





38 


1.552 ±0 


Q41+0.12 


1.80 ± 2.87l|^^ 


DIS 





50- 


-0 


58 





54 


1.238 ±0 


Q2q+o.io 


4.80 ± 0.441°:^? 







58- 


-0 


66 





62 


1.298 ±0 




5.53 ± QMtoll 







66- 


-0 


74 





70 


1.324 ±0 


Q24+0.09 


6.51 ± 0.22l°:^^ 







74- 


-0 


82 





78 


1.276 ±0 




7.65±0.19l°:^8^ 







82- 


-0 


90 





86 


0.975 ± 


0191°:°^ 


8.40 ± 0.191°:^^ 







90- 


-1 


00 





93 


0.248 ± 


oo8l°:J° 


7.99 ± 0.28l^:^^ 


high-Q^ 





20- 


-0 


50 





38 


1.786 ±0 


Q7C+0.13 
"-""-0.15 


-0.53 ± 4.56l|^^ 


DIS 





50- 


■0 


58 





54 


1.449 ± 


0561°:°^ 


5.79 ± 0.741°:?^ 







58- 


-0 


66 





62 


1.309 ±0 


0451°:!^ 


4.66 ± 0.481°:^^ 







66- 


-0 


74 





70 


1.374 ±0 


0441°:^^ 


6.76 ± 0.401^:1 







74- 


-0 


82 





78 


1.283 ±0 


0391°:°? 


7.02 ± 0.33l|]:^^ 







82- 


-0 


90 





86 


0.915 ±0 


nqo+O.OS 
"-''-'^-0.09 


8.29 ± 0.36l[J:?} 







90- 


-1 


00 





93 


0.217 ±0 


Qi 0+0.10 
"-'-'-"-'-0.08 


7.93 ± 0.53l°:| 



Table 6: The intercepts a and slopes h from the exponential parameterization of 



the differential cross defined in Section for the photoproduction and three DIS 
samples. Statistical uncertainties are listed first, followed by systematic uncertain- 
ties, not including an overall normalization uncertainty of 2.1% on the intercepts 
in DIS, nor an additional uncorrelated uncertainty of 5.1% on the photoproduction 
intercepts. The systematic uncertainties are largely correlated between all points. 
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e' 




Figure 1: a) HERA ep scattering event with the final-state baryon in the proton- 
fragmentation system, X . b) Leading baryon production via an exchange process. 
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Figure 2: Side view of the proton beamline downstream from the ZEUS inter- 
action region. The protons are moving from right to left. The labels for the HERA 
components, e.g. B^l, indicate the horizontal distance in meters from the interac- 
tion point. The horizontal and vertical axes are not to scale. The proton beam is 
bent upward by approximately 6 mrad by the dipole magnets B67-B77 near Z=+70 
m. The FNC is located on the zero-degree line at Z=-hl05.5 m. S1-S6 indicate the 
locations of the ZEUS leading-proton spectrometer stations [26]. 
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c 




— 1 


B — 









Figure 3: Diagram of the FNC/FNT assembly. The thick horizontal lines show 
the 5 cm vertical segmentation of the front part of the FNC. The hole through 
the third and fourth towers from the top allows the proton beam to pass through 
the calorimeter. The 17 x 15 grid of small squares shows the fingers of the FNT 
hodoscopes. The irregular curve shows the geometric aperture defined by upstream 
beamline elements, and the bullet (•) shows the zero-degree point. 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 4: The kinematic regions in (a) p^ and (b) t covered by the angular 
acceptance of the FNC (On < 0.75 mrad) are shown as shaded bands. 
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4 6 8 10 12 14 

X(cm) 

Figure 5: Scatter plot of reconstructed positions in the FNT. The irregular curve 
is the geometric aperture defined by upstream beamline elements as modeled in the 
Monte Carlo. The zero-degree point is at X = 12.5 cm, Y = 8.3 cm. The deficit of 
events observed near Y ~ 10.5 cm is due to dead material in the LPS. 
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Figure 6: Vertical slice of the data in Fig. 0/or 7 < X < 10.5 cm. The smooth 
curve is the generated Monte Carlo distribution; the histogram is the Monte Carlo 
distribution of reconstructed events after the dead material simulation. The peak of 
the smooth curve, determined from a fit to the data, corresponds to the zero-degree 
point in the vertical plane. 
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Figure 7: The distribution {1/ <Ji-ac)d<^iM/ dxi for neutrons in DIS with scattering 
angles On < 0.75 mrad, corresponding to the kinematic range p\. < 0.476 x| GeV"^ . 
The statistical uncertainties are smaller than the plotted symbols; the shaded band 
shows the systematic uncertainties. The band does not include the overall normal- 
ization uncertainty of 2. Wo. 
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Figure 8: The p\. distributions for DIS. Note the logarithmic vertical scale 
and the varying p"^ ranges. The statistical uncertainties are shown by a vertical 
error bar; in most cases they are smaller than the plotted symbol. The systematic 
uncertainties are not shown. The line on each plot is the result of a fit to the form 
da^^/dpip oc exp(— 
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Figure 9: The intercepts a = {1/ ainc)(Po-LN / dxLdp'^\p2^^Q versus Xl from fits 

of the distributions in DIS to the form dau^/dp\. = aexp(— over the p^ 
ranges shown in Fig. O The statistical uncertainties are shown by a vertical error 
bar; in most cases they are smaller than the plotted symbol. The shaded band shows 
the systematic uncertainties. The band does not include the overall normalization 
uncertainty of 2.1%. 
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Figure 10: The exponential slopes b versus Xl from fits of the p^ distributions in 
DIS to the form dau^/dp^ = aexp(— over the p^ ranges shown in Fig. 0. The 
error bars show the statistical uncertainties; the shaded band shows the systematic 
uncertainties. 
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Figure 11: The ratio a/b versus xl from fits of the distributions in DIS 
to the form dai^^/dp^ = aexp(— over the ranges shown in Fig. 0. The 
error bars show the statistical uncertainties; the shaded band shows the systematic 
uncertainties. The band does not include the overall normalization uncertainty of 
2.1%. 
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Figure 12: The xl distributions for the photoproduction and three DIS suhsam- 
ples. The error bars show the statistical uncertainties. The common systematic 
uncertainties are shown as a shaded hand. There is an overall normalization un- 
certainty of 2.1% for the DIS data, and an additional uncorrelated uncertainty of 
5.1% for the photoproduction data. 
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Figure 13: The exponential slopes b for the photoproduction and three DIS 
subsamples. The points are slightly offset horizontally for clarity. The error bars 
show the statistical uncertainties. The common systematic uncertainties are shown 
as a shaded band. 
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Figure 14: The ratio of the Xl distributions for photoproduction and DIS, nor- 
malized as described in the text, for neutrons with scattering angles On < 0.75 mrad, 
corresponding to the kinematic range p\ < 0.476 x| GeV"^ . The error bars show the 
statistical uncertainties; the shaded band shows the systematic uncertainties. The 
band does not include the overall normalization uncertainty of 5.1% on the ratio. 
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Figure 15: The p\ distributions for photoproduction and DIS. They are each 
normalized to unity at p\ = 0. Note the logarithmic vertical scale and the vary- 
ing p\. ranges. The statistical uncertainties are shown by a vertical error bar; 
in some cases they are smaller than the plotted symbol. The systematic uncer- 
tainties are not shown. The lines on each plot are the results of fits to the form 
(icTLN / dp"^ oc exp {—bp"^) . 
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Figure 16: The differences between the exponential slopes in photoproduction and 
DIS, Ab, versus xl. The error bars show the statistical uncertainties; the shaded 
band shows the systematic uncertainties. 
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Figure 17: The normalized differential distributions / (y\nc) d<^h^/ d^L for leading 
protons and neutrons in the range < 0.04 GeV^. For the leading-neutron points, 
the statistical uncertainties are shown by a vertical error bar; in most cases they 
are smaller than the plotted symbol. The shaded band shows the systematic uncer- 
tainties; the band does not include the overall normalization uncertainty of 2.1%. 
For the leading-proton points, the inner error bars show the statistical uncertainty 
only; the full error bars show the statistical and systematic uncertainties added in 
quadrature. 
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Figure 18: The exponential slopes of the p'^ distributions for leading neutrons and 
protons in DIS. For the leading-neutron points, the error bars show the statistical 
uncertainties; the shaded hand shows the systematic uncertainties. For the leading- 
proton points, the inner error bars show the statistical uncertainty only; the full 
error bars show the statistical and systematic uncertainties added in quadrature. 
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Figure 19: Comparison to Monte Carlo models of DIS: (a) xl distributions, (b) 
intercepts, and (c) exponential slopes. The error bars show the statistical uncertain- 
ties; the shaded bands show the systematic uncertainties. The bands do not include 
the overall normalization uncertainty of 2.1% for the Xl distribution and intercepts. 
The curves are from the DIS Monte Carlo models Rapgap [24] and Lepto [28]. 
The curves labeled RAPGAP std. frag, and LEPTO std. frag, are the models incor- 
porating only standard fragmentation. The curve labeled RAPGAP-vr includes also 
Pomeron and pion exchange; the curve labeled LEPTO-SCI is the model including 
soft color interactions. 
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Figure 20: The measured exponential slopes b compared to the predictions of mod- 
els of one-pion exchange alone. The error bars show the statistical uncertainties; 
the shaded band shows the systematic uncertainties. The curves show predictions 
from the models discussed in the text [6-8]. 
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Figure 21: Ratio of photoproduction and DIS Xl distributions. The error bars 
show the statistical uncertainties; the shaded band shows the systematic uncertain- 
ties. The band does not include the overall normalization uncertainty of 5.1% on the 
ratio. The dashed curve shows the neutron loss from a model of rescattering [I4]; 
the solid curve is this model corrected for the different W dependence of the pion 
cross section in DIS and photoproduction. The dot-dashed curve shows the neutron 
loss from another absorption model [13]; the dotted curve is this model corrected for 
the different W dependence of the pion cross section in DIS and photoproduction. 
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Figure 22: Xl distribution for photoproduction. The error bars show the statistical 
uncertainties; the shaded band shows the systematic uncertainties. The band does 
not include the overall normalization uncertainty of 5.5%. The dashed curve shows 
the prediction of a model with enhanced neutron absorption and migration for pion 
exchange only [15]. The solid curve shows the same model including also p and 02 
exchanges [16]. 
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Figure 23: (a) Exponential slopes b for DIS, and (b) difference of exponen- 
tial slopes b for photoproduction and DIS. In both plots the error bars show the 
statistical uncertainties; the shaded band shows the systematic uncertainties. The 
dashed curves are from a model [16] based on pion exchange with enhanced neutron 
absorption and migration; the solid curves include also p and 02 exchanges. 
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